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Abstract 



The light quark masses are determined using a new QCD Finite Energy Sum Rule (FESR) in the pseu- 
doscalar channel. This FESR involves an integration kernel designed to reduce considerably the contribu- 
tion of the (unmeasured) hadronic resonance spectral functions. The QCD sector of the FESR includes 
perturbative QCD (PQCD) to five loop order, and the leading non-perturbative terms. In the hadronic sec- 
tor the dominant contribution is from the pseudoscalar meson pole. Using Contour Improved Perturbation 
Theory (CIPT) the resuhs for the quark masses at a scale of 2 GeV are m„(Q = 2 GeV) ^ 2.9 ± 0.2 MeV, 
TOd(Q = 2 GeV) = 5.3 ± 0.4 MeV, and m,(Q = 2 GeV) = 102 ± 8 MeV, for A = 381 ± 16 MeV, corre- 
sponding to Q!s(M^) = 0.344 ± 0.009. In this framework the systematic uncertainty in the quark masses 
from the unmeasured hadronic resonance spectral function amounts to less than 2 - 3 %. The remain- 
ing uncertainties above arise from those in A, the unknown six-loop PQCD contribution, and the gluon 
condensate, which are all potentially subject to improvement. 
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The QCD running quark masses, together with the running couphng constant, are the most important 
parameters of the strong sector of the Standard Model of Particle Physics. Many attempts have been made 
over the years to extract the values of the quark masses from hadronic observables in various theoretical 
frameworks. This pursuit started with current algebra, and later chiral perturbation theory, which fixes 
light quark mass ratios ,1 . This was followed by QCD sum rules Il]-[5], and lattice QCD [3]. In the 
framework of QCD sum rules the ideal Green's function is the pseudoscalar correlator, as it involves the 
quark masses as an overall multiplicative factor (with subdominant and small quark mass corrections), 
and the well known pseudoscalar meson pole. Unfortunately, the hadronic resonance spectral function is 
not directly known from experiment. In fact, except at best for the masses and widths of the first two 
radial excitations in the non-strange and strange sector, nothing is known about the actual shape of these 
spectral functions. This information is hardly enough to realistically reconstruct them; inelasticity and 
non-resonant background would remain unknown. We describe here a recent attempt [5]- [6] to reduce 
considerably this systematic uncertainty by incorporating a specific analytic kernel in Cauchy's theorem 
used to derive the FESR. This kernel takes the form of a second degree polynomial which vanishes at the 
peaks of the two radial excitations of the ground state pseudoscalar meson (pion or kaon). The idea of 
including polynomial kernels in FESR is not new; several different forms have been used in the past for 
various reasons and in various channels. What is new, though, is the requirement that the kernel vanishes 
at the peaks of the resonances entering the experimentally unknown spectral functions. This kernel does 
distort whatever model one uses as a guess for the resonance spectral function. But then, it also distorts 
the pseudoscalar meson pole, and the QCD contribution. Since Cauchy's theorem is still valid, so is the 
distorted FESR. Remarkably, this procedure not only reduces considerably the importance of the unknown 
sector, but at the same time it leads to an exceptionally broad region of stability against changes in sq, 
the upper limit of integration in the FESR, which in CIPT is Sq — 1 ~ 4 GeV'^. We consider the correlator 
of (light) axial-vector current divergences 



Uq')\i = ^ (m, + m,) / d\ e'^^x < T\{Mx),jlm\ > \i , (1) 



where j5{x)\j —: qj{x) 475 qi{x) : , and i,j are quark flavours. To simplify the notation we shall use in the 
sequel rrij + mi = m. We define the integration kernel 

^5(5) = 1 - ao s - ai , (2) 

where qq, and ai are free parameters to be fixed by the requirement that A5(Mf ) — A5(M|) = 0, with 
Mi_2 the masses of the two resonances in the pseudoscalar channel, 7r(1460) and 7r(1830) for the non- 
strange channel, and K(1460) and K(1830) for the strange channel. Invoking Cauchy's theorem for the 
second derivative of ^'5(9^) one finds [6] 



where 

/ ip- \ ^2 

F(s) = -s ( So - ao y - ai y j + y - ao - ai ■ (4) 

The function "05 '^'^^(9^) is currently known up to five-loop order in PQCD [7j. We define the left hand 
side of Eq.(3) as 
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Figure 1: Strange quark mass at 2 GeV for Aqcd = 365 (397) MeV, curves (a) and (b), respectively. 



S5{so)\qCD = "TT- f ^^"^5 {s)\qCD [F{s) - F{so)], 



c(U\ 



(5) 



and obtain in PQCD and in the framework of CIPT [6] 



S5{so)\pQCD 



F{x)-F{so) [as{x)Yexp 



M=0 



2i J2 IM / dx'[as{xX'^^ 



(6) 



where the RGB for the mass and couphng have been used, as{x) = as{x)/Ti, and where 



F{x) = Y.^-rb^s 



N iNx 



(7) 



iV=l 



with the angle x defined through s = sq e*^, and x e (— tt, tt). 



The constants above are {np — 3): — Coi, Ki = Cn + 2Ci2, K2 = C21 + 2C22, ~ C31 + 2C32, 
= C41 + 2C42, with Coi = 6, Cii = 34, C12 = -6, C21 = -105 C(3) + 9631/24, C22 = -95, 
C23 = 17/2, C31 4748953/864 - 7rV6- 91519 C(3)/36 + 715 C(5)/2, C32 = -6 [4781/18- 475 C(3)/8], 
C33 = 229, C34 = -221/16, C41 = 33532.26, C42 = -15230.6451, C43 = 3962.45493, C44 = -534.052083, 
C45 = 24.1718750, and C,{x) is Riemann's zeta function. Finally, hi = — (sq — aoSg/2 — aiSQ/3), 62 = 1/2, 
63 = — ao/6, and 64 = — ai/12. Regarding the value of Aqcd entering Q!s(so), it can be extracted 
from the strong coupling obtained from r-decay . We use the latest high precision result of [S], i.e. 
as{M^) = 0.344 ± 0.009, leading to Aqcd = 365 - 397 MeV. 



The light-quark condensate contribution to the left hand side of Eq. (3) is 
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m^(.So) _ , . 
d5[So)\<qq> = -2 2 — \^sqq}\^o 
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dx e 



F(x)-F(so) 
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1 + Y"-s(x) 



X exp 
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Figure 2: Up quark mass at 2 GeV for Aqcd — 365 (397) MeV, upper (lower) curve, respectively. 



and the contribution of the gluon condensate becomes 



^5(so) 



<G2> 



1 m^iso) a. 



-G )Uo - 



dx e 



16 , , 121 , , 
F{x) - F{so) 1 + yas(^o) + — as(x) 



X exp 



2iJ2jm f dx'[as{xX'^' 



(9) 



where the scale /ig ~ 1 GeV^ appears in connection with the removal of logarithmic quark mass singularities 
(see [3])- The dimension six four-quark condensate, and the higher order quark mass corrections turn out 
to be negligible in this application. We shall use (^G^) ~ 0.06 GeV^ and (qq) ~ (-250 MeV)^ 0. CIPT 
[To] has been shown to provide better convergence than Fixed OPT in the QCD analysis of the vector and 
axial- vector correlators in tau-lepton decay [9]- [10]. We found this to be also the case in the determination 
of the light quark masses [6] . Unlike the case of FOPT, where as(so) is frozen in Cauchy's contour integral 
and the RG is implemented after integration, in CIPT as and the quark mass are running and the RG is 
used before integrating. This is done through a single-step numerical contour integration and using as input 
the strong coupling and the quark mass obtained by solving numerically the RG equation. Notice that the 
invariant quark mass enters the expression of the running quark mass as an overall multiplicative factor; 
this invariant mass is what one determines from the FESR. Turning to the hadronic resonance sector, we 
have used a parametrization involving two Breit-Wigner forms normalized at threshold according to chiral 
perturbation theory, as first proposed in [llj . It is unnecessary to improve on this parametrization as the 
presence of the integration kernel in the FESR reduces this contribution to the level of 2 - 3 %. In Figs. 
(l)-(3) we show the results for the strange, the up, and the down quark masses at Q = 2 GeV. The region 
of stability, sq — 1 — 4GeV^, is exceptionally wide, with a remarkably low spread of 4 - 5 % for the up 
and down quark masses, and 1 - 2 % for the strange quark mass. To achieve a reasonable error we have 
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Figure 3: Down quark mass at 2 GeV for Aqcd — 365 (397) MeV, upper (lower) curve, respectively. 



allowed for a ±30 % uncertainty in the hadronic resonance parametrizations, a factor two uncertainty in 
the value of the gluon condensate, and have assumed the (unknown) six-loop contribution to be equal to 
the five-loop one. In order to disentangle the strange quark mass from the sum (m^ +mq) (q=u,d) we have 
used the chiral perturbation theory ratio rns/niud = 24.4 ± 1.5, with niud = (wu + 'md)/2. Similarly, 
to disentangle e.g. the up quark mass we used niu/md = 0.553 ± 0.043. This leads to the values 



m„(2 GeV) = 2.9 ± 0.2 MeV , 



(10) 



TOd(2GeV) =5.3 ± 0.4 MeV 



m,(2GeV) = 102 ± 8 MeV 



(11) 
(12) 



m„ -I- rud 



= 4.1 ± 0.2 MeV 



(13) 



rrid + rriu 
rud - m„ 



0.29 ± 0.05 



(14) 



These results are in agreement with some recent lattice determinations [4]. 
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